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ABSTRACT 

We present a study of the kinetic energy dissipation in interstellar cloud collisions. The 
main aim is to understand the dependence of the elasticity (defined as the ratio of the final to 
the initial kinetic energy of the clouds) on the velocity and mass ratio of the colliding clouds, 
magnetic field strength, and gas metallicity for head-on collisions. The problem has been 
studied both analytically and via numerical simulations. We have derived handy analytical 
relationships that well approximate the analogous numerical results. The main findings of this 
work are: (i) the kinetic energy dissipation in cloud collisions is minimum (i.e. the collision 
elasticity is maximum) for a cloud relative velocity v r ~ 30 km s _1 ; (ii) the above minimum 
value is proportional ZL^ where Z is the metallicity and L c is the cloud size: the larger is 
Zl? c the more dissipative (i.e. inelastic) the collision will be; (Hi) in general, we find that the 
energy dissipation decreases when the magnetic field strength, and mass ratio of the clouds are 
increased and the metallicity is decreased, respectively. We briefly discuss the relevance of this 
study to the global structure of the interstellar medium and to galaxy formation and evolution. 



Subject headings: Hydrodynamics - Shock waves - ISM: kinematics and dynamics - Numerical 
methods 



Introduction 



The complex morphology and the variety of scales observed witness to the turbulent nature of the 
ISM. Nevertheless, the vast majority of the dynamic models of the ISM (but see Norman & Ferrara 1996, 
Vazqucz-Scmadcni et al. 1995) describe this medium in terms of dense, neutral clouds (CNM) moving into 
a smooth, diffuse intercloud medium (WNM) by which are pressure confined. If we assume that the cloud 
velocities are randomly distributed, a cloud undergoes a supersonic collision approximately every 10 6 — 10 7 
years. Therefore, in the Galaxy, assuming a volume of the gaseous disk to be ~ 78 kpc 3 , about 0.3 — 3 cloud 
collisions occur every 100 years (we have assumed that a line of sight through the Galactic disk intersects 
~ 7 clouds per kpc [Spitzer, 1978]). 

Cloud collisions play an important role in the evolution of a galaxy. For example, the large gas 
compression following an inelastic collision might be expected to enhance the star-formation rate (Gildcn, 
1984). Inelastic cloud collisions are important for the global energy budget of the ISM: as the energy input 
from SNe and HII regions accelerate the interstellar gas, the cloud kinetic energy is dissipated by cloud 
collisions (Spitzer, 1978). Also, since the vertical distribution of the HI in the Galaxy depends on the 
turbulent motions of the clouds (McKee 1990, Ferrara 1993), the structure of the gas is strongly regulated 
by the dissipation of bulk motions. 

In addition, the buildup of the mass spectra of an ensemble of clouds is partially affected by cloud 
collisions (Oort 1954, Field & Saslaw 1965, Field & Hutchings 1968, Penston et al. 1969, Cowie 1980, 
Hausman 1982, Pumphrey & Scalo 1983, Struck-Marccll & Scalo 1984, Fleck 1996). The basic idea is that 
small clouds, due to inelastic collisions, coalesce to form larger ones which eventually collapse to form 
stars. Klein, McKee & Woods (1995) have recently pointed out that the outcome of an inelastic collision 
may be shattering rather than coalescence due to the fast growth of a nonlinear instability, the so-called 
bending mode instability. This result has a significant impact on the cloud buildup models. Moreover, the 
morphology of the clouds after the collision event is strongly dependent on both the elasticity (defined as 
the ratio of the final to the initial kinetic energy of the clouds) and the presence of a magnetic field. 

The rationale for the assumption of inelastic collisions and subsequent coalescence of the colliding 
clouds commonly adopted by the buildup models can be summarized as follows. Interstellar clouds motions 
are highly supersonic; therefore, during collisional events, two shock waves arise and propagate from the 
contact discontinuity (CD) heating the gas; this hot gas, being optically thin, tends to cool quickly due to 
radiative losses. For a standard diffuse cloud, the typical cooling time of the shocked gas is much less that 
the characteristic timescale of the collision. This justifies the inelastic approximation. Nevertheless, this 
assumption might be rather crude for small clouds and/or primordial galaxies where the gas metallicity is 
low and therefore the cooling time is longer. 

In this work we will investigate the physics of interstellar cloud collisions, with particular emphasis 
devoted to the dependence of the kinetic energy dissipation on the parameters of the problem, such as 
the cloud relative velocity, gas metallicity, cloud mass ratio, and magnetic field strength. The final goal 
is to identify the regions of the parameter space in which the collisions are elastic (adiabatic), inelastic 
(isothermal) or radiative. 

The plan of the paper is as follows. In § II we discuss the main heating/cooling processes of the ISM 
and calculate its thermal balance. In § III, we will present a first analytical approach that will allow us 
to find approximate relationships between the elasticity and the cloud parameters (relative velocity, sizes, 
gas metallicity, dust-to-gas ratio). These analytical results should be useful for future works concerned 



with a statistical approach to the dynamics of the galactic ISM, the mass spectrum of diffuse clouds, and 
star formation processes. In § IV, the results of a series of numerical simulations of head-on collisions are 
considered; in § V we will discuss the results and the dependence of the elasticity on the various parameters 
of the collision comparing the analytical and numerical approaches. A brief summary closes the paper. 



2. Thermal balance of ISM 



In this Section we calculate the thermal equilibrium of the diffuse interstellar medium that is illuminated 
by the local interstellar far-ultraviolet and X-ray radiation field and permeated by the cosmic-ray flow. We 
follow closely the work of Wolfire et al. (1995), which have incorporated the photoelectric heating from small 
grains and PAHs (with a distribution of sizes in the range 3-100 A). Considering a realistic cooling/heating 
function is important for two reasons: (i) we want to study the energy dissipation in collisions occurring in 
a two-phase medium; (ii) the elasticity of the collision is closely connected to the emission processes. In 
addition we would like to explore the effects of gas metallicity, Z, and dust-to-gas ratio, D/G (normalized 
to the local value), on the equilibrium states of the ISM. A summary of the dominant cooling processes 
included in our function is given in Table 1. The Table is intended to complement the analogous one in 
Wolfire et al. (1995). 

In order to obtain an analytical expression for the fractional ionization of H and He, and the electron 
density n e , we have introduced the approximations described below. The motivation for this choice is to 
decrease the computational time of the numerical simulations which make use of the cooling function. 
We assume the ionization equilibrium and we neglect: recombination on dust grains, He dielectronic 
recombination, ionization due to hydrogen collisions, secondary ionization by cosmic-rays; in addition, 
we assume x^+ = n H +/n = Xtf e + = n He + /n{He), where n = nn + n H +, n(He) = nn e + n He + and 
n(He)/n = 0.1. With these approximations the equation of the ionization equilibrium can be solved 
analytically to obtain an expression for n e (n,T). 

In Figure we show the phase diagram (pressure versus density) of the gas for different values of the 
metallicity and of the dust-to-gas ratio. When Z (which we will assume to be equal to D/G throughout the 
paper) is decreased the pressure range in which a multi-phase medium can exist becomes wider; in addition 
the mean equilibrium pressure and density of the clouds increase. The photoelectric heating of PAHs and 
small grains dominate on the cosmic-ray and X-ray heating both in the intercloud medium (WNM) and in 
the clouds (CNM). The WNM is cooled mainly by the emission of Lya, CII (158 /im) and 01 (63 /im). The 
CNM is cooled mainly by CII (158 /im). 



3. Analytical approach 



Before we describe the results of the numerical simulations, we present here a simple analytical model 
for the collision. As already stated in the Introduction, the motivation is twofold: on the one hand we can 
easily isolate the most important physical processes; on the other hand simple approximated relationships 
can be found that can be used to explore a wide region of the parameter space. 



As a first approach to a parametric study of the energy dissipation in interstellar cloud collisions we 
focus on a ID model, which allows us to consider centered collisions only; the study of off-center collisions is 
devoted to a future paper. We will initially assume that the two clouds are identical and we will release this 
assumption in Sec. 3.4. In this simple scenario, clouds are modeled as cylinders of base area S c and initial 
length L c . The time origin t = is defined as the moment at which the bases of the cylinders coincide in a 
plane, i.e. the contact discontinuity. We use a reference frame (LSR) with the ir-axis along the symmetry 
axis of the cylinders and the y — z plane coinciding with the CD. In the LSR, v c is the cloud velocity, Vf 
is the shock velocity (see below), and vcd the velocity of the CD (for the collision of two identical clouds, 
vcd — 0, and for the symmetry of the problem we can consider the evolution of one cloud only); the clouds 
have initial density no, and a mass Mq = S c L c pmHno, where p is the mean molecular weight of the gas and 
ran is the proton mass. 

Three main evolutionary stages of the collision can be identified: (Stone, 1970a, b): (1) a compression 
phase: two shocks propagate from the CD into the clouds converting the kinetic energy of the unperturbed 
gas into internal energy of the shocked gas. At the end of this phase the clouds are crushed into a disk 
whose thickness depends on the entity of radiative losses. (2) an expansion phase: when the shock reaches 
the rear side the cloud expands because of the decrease of the ram pressure of the incoming flow. The shock 
will propagate into the intercloud medium and a rarefaction wave will propagate backwards in the cloud. 
(3) a collapse phase: the rarefaction wave brings the pressure of the cloud below the ambient pressure: 
when the rarefaction wave is reflected by the cloud leading edge it becomes a compression wave thus halting 
the expansion. In this phase the cloud rear edge becomes Raylcigh- Taylor unstable. 

Depending on the values of N coo i — noV c t C ooU the column density of the post-shock radiative region 
(where t coo i is the typical cooling time of the shocked gas), and N c = uqL Ci the total cloud column density, 
the collision will be defined to be: (a) adiabatic: N coo i 3> uqL c (t CO ol 3> t co u); (b) isothermal: N coo i <C n§L c 
(tcooi < tcou); (c) radiative: N coo i ~ n L c (t coo i ~ t cM ); where t coU = L c /v c is the typical collision time. In 
the following we will briefly discuss the limiting cases (a) and (b) and present in more detail the radiative 
case. 



3.1. Adiabatic collisions 



The shock wave generated by the collision event will propagate from the CD through the cloud with 
a velocity Vf in the LSR and v s = v c + Vf with respect to the unperturbed gas; the shocked gas is at rest 
in the LSR and moves with velocity vj with respect to the shock front. For a strong shock in a perfect 
gas with 7 = 5/3, the density pi, the pressure, p\, the temperature, T\ and the sound speed, Ci, of the 
postshock gas are: 

^ = ^4, %*\m\ ^A^^f.. (1) 

Po Vf Pq 4 T 16 4 

At the end of the compression phase, t ~ t C oll, the gas density is enhanced by a factor 4, the shock velocity 
is v s — (4/3)v c , and the cloud length is L ~ L c /A. By definition, for an adiabatic collision the radiated 
energy, E r , is zero: 

— = 0, (2) 



where El = Mv 2 is the initial kinetic energy of the clouds. The elasticity of the collision is: 



e -Ei~Avl- 1 - (3) 



where Ej, = Mv\ ~ (9/4)MC 2 is the kinetic energy of the cloud in the expansion phase. 



3.2. Isothermal collisions 

In an analogous manner, for isothermal collisions we find: 

^? = ?^» £ = £-!■ (4) 

■Po Pa v f T 6 

As a consequence, when the shock reaches the rear side of the cloud, the cloud size is reduced to L ~ L c /A4 2 
and the shock velocity is v s — [A4 2 /(A4 2 — l)]v c ~ v c . In this case the initial kinetic energy will be almost 
completely dissipated and e ~ Mr 2 ~ 0. 



3.3. Radiative collisions between identical clouds 

In the intermediate cases in which N coo i ~ no£ c (radiative collisions), it can be shown that the 
elasticity of the collision depends on a single parameter r\ ex N coo i/(noL c ); we will refer to 77 as the "elasticity 
parameter" . After the shock has reached the rear side of the cloud, the cloud begins an expansion phase 
that we will assume to be adiabatic. This assumption is reasonable since the typical expansion timescale is 
much shorter than t CO ol- 

At this time there are two physically different regions in the cloud interior: (i) an adiabatic layer 
behind the shock front; we will assume that the column density of this region is proportional to the column 
density N coo i of the radiative region of the shock, (ii) an isothermal region in which the temperature of the 
gas is approximately in thermal equilibrium. 

As we have seen in the previous Section, the major contribution to the radiated energy comes from the 
isothermal region, whereas the kinetic energy is essentially stored into the adiabatic region. If we define the 
"elasticity parameter" of the collision , 77, as 

I* cool / r \ 

r] = a — , (5) 

n L c 

where a is defined by N ac i = a ■ N coo i (i.e. we assume that the column density of the adiabatic region is 
proportional to the cooling column density) we obtain the following relationship for the radiated energy 
and collision elasticity e: 

l-e=!r - (!-»?). ( 6 ) 



El 

4 



m (7) 



obviously, 77 <§C 1 corresponds to isothermal collisions; 77 = 1 denotes adiabatic collisions 
From a dimensional analysis, we find: 

2 l +e 'C{x e e q ,n-1 7 T) 



t coo i = -(1.1 + x e e q ) f ,^ eq _ M ^ . (8) 



where ST = T - T e i. Substitution of eq. (|) into eq. (§) yields, 

v s n 5T 

71 ""W (9) 

Behind the shock front it is |£| = |?1qA(T) — nor(T)| ~ nf ) K(T) 1 because the radiative losses dominate 
over the heating and ST oc v^ for the adiabatic Rankine-Hugoniot jump conditions. It follows that 

3 

VK—j-£j sv (10) 

noL c A(v s ,Z) 

In Figs. and we show (solid lines) the dependence of (1 — 77) = E r /E l k on the relative velocity, v r , of 
the clouds for different values of L c , and Z as given by eq. [Tol In Sec. 5 we will show how the above 
analytical solutions approximate quite well the analogous results obtained from the numerical simulations. 
We have assumed that v s ~ v c = v r /2 (this hypothesis holds exactly for a completely inelastic collision, see 
Sec. 3.2). The number density no and temperature To of the unperturbed clouds are computed imposing 
pressure equilibrium with the intercloud medium; if D/G = Z = 1 we assume P/k — 2300 cm~ 3 K for the 
equilibrium pressure. If Z 7^ 1, the equilibrium pressure is computed fixing the CNM temperature to the 
value found with D/G = Z = 1, To ~ 48 K; thus, if the metallicity decreases, clouds become more dense 
and the equilibrium pressure increases. The motivation for this choice is that the equilibrium temperature 
of the clouds is poorly sensitive to the value of the metallicity and pressure of the ISM. 



3.4. Radiative collisions between different clouds 



We now consider colliding clouds with different velocities and sizes; the indexes 1 and 2 refer to 
parameters of the left and right cloud, respectively. We define the nondimcnsional parameters: 



icl = Ml g 



(11) 



Since we have assumed the same number density no for the two clouds the size ratio is the same as the 
mass ratio; vcd, is computed imposing the same ram pressure on the two sides of the discontinuity. Until 
the shock wave is inside the smaller cloud we have: 



\vcd\ 



7r(0-l). (12) 



2 
In the frame of reference comoving with the CD (DSR), the cloud velocity will be 



y = M(/? + l), (13) 



where v r is the relative velocity of the clouds. This phase of the collision (in the DSR) is the same as for 
two identical clouds. The kinetic energy in DSR, E' k , is related to the kinetic energy in the LSR, E k , by the 
relationship: 

El = f(K,0)-Bt, (14) 



where 



/(«,/?) 



/, =f(n,P)-El 

(0+1) 2 (k + 1) 



(15) 



4(1 + k/? 2 ) ' 

obviously if k = = 1 then /(l, 1) = 1. Since E r is the same in both reference frames, then the elasticity 
in the LSR is just the product of the elasticity in the DSR and /(«, 0). 

As the shock wave exits from the smaller cloud, the discontinuity starts to move (in the DSR) under 
the effect of the ram pressure of the bigger cloud. The dynamics becomes rather complex since vqd depends 
on the ram pressure (oc v 2 ) and v s depends on vcd as v s (t) = v s (0) — vco{t) (in the DSR). During this 
phase the shock velocity is not constant and the radiated energy must be derived by integration: 



E r 

El 



1 



S c N c v 2 c 



dt n n AS r L r = 



N c - aN cool 

N~ r 



1 



N c -aN cool 



Ncv 2 c j 



S(l-Tl), 



dN vi = 



where 



N c -aN cool 



dN v 2 = £,v 



(16) 



(17) 



N c - aN coo i J a 

where L r is the length of the post-shock radiative region, aN coo i is the column density of the adiabatic 
region at the moment in which the shock exits from the cloud. If v s (t) = cost., then £ = 1 and we recover 
the case of identical clouds: E r jE\ = (1 — 77). After some algebra we find (in the DSR): 



E r 

El 






where 



Vi 



N C ool(Vs 



m = 



N C ool(Vs2) 



(18) 



(19) 



N cl ' " N c2 ' 

v s i is the shock velocity exiting from cloud i, N C i is the column density of cloud i. If we assume that 
vqd ~ w cm/2, where vqm is the velocity of the center of mass, at the end of the collision, we have: w s2 ~ v c 
and v s i ~ v c (k + 3)/[2(k + 1)]. We derive ^ from an energy conservation argument: if all the kinetic energy 
relative to the center of mass is dissipated we obtain 771 = 772 = 0. In this case E r jE k = [2 + £(k — 1)]/(1 + k) 
is equal to 1 — (Eq m I 'E k ) = 4k/(k + l) 2 , where Ecu is the kinetic energy of the center of mass. This 
allows us to conclude that £ = 2/(k + 1); finally in the LSR we have: 



l-e=§=/(«,/3) 



2-£ 



(l-»») 



< 



(1-m) 



(20) 



3.5. Lateral outflow 



In a realistic situation clouds have a finite size perpendicular to the collision direction (x-axis); thus 
the gas will be free to expand in the directions parallel to the shock front, possibly generating a lateral 
outflow. We will compute the amount of mass lost in the lateral outflow analytically. To this purpose, we 
integrate the flux of mass from the cloud surface on the typical collision time t co u — L c /v c . The typical 
expansion velocity of a gas in the vacuum is proportional to the sound speed; as a consequence, most of the 
mass loss comes from the hotter post-shock adiabatic region. During the compression phase the amount of 
mass ejected in the lateral outflow will be: 

Mout = ^L y L ad v out t coU Ap , (21) 

where L y is the cloud diameter (clouds are modeled as cylinders), L ad if the length of the adiabatic region, 
Vcmt is the outflow velocity, and 4po is the post-shock gas density. After some algebra, and recalling that 
v ut — v s and Ap L ad = aN cooh we find: 

*£"*■ (22) 

where rjy = N coo i/(ri()Ly). The kinetic energy, Ek.out ck M out v^ ut1 lost in the lateral outflow is: 

-E^k.out 



E' k 



Vy (23) 



The kinetic energy lost in the lateral outflow depends both on the elasticity and the geometry of the 
collision. Since r\ y = rj(L c /L y ), (a) for spherical clouds (L c /L y = 1), the kinetic energy lost in the lateral 
outflow is proportional to the elasticity parameter (i.e. is zero in a perfectly inelastic collision and 1 in a 
perfectly elastic one); (b) for a fixed of r/, the energy loss is smaller for flattened clouds (L y ^> L c ). 



We can rewrite eqs. (9) and (23) including the energy losses in the lateral outflow 



1 


^k 


which is valid for a head-on collision between identical clouds, 


l-e=f^/(«,/3) 


— i l - Vy i - m + —5- 

1 + K 1 + K 



(24) 



-(1 - Vy )(l - m ) 



(25) 



for a head-on collision between different clouds. These results are obtained supposing that the lateral 
outflow is adiabatic. 



Numerical results 



In this Section we study cloud collisions via numerical simulations adopting the same ID model as in 
the analytical approach. The numerical code is based on a shock-capturing scheme (Yee, 1989) suitable 
to resolve the hypersonic shocks that arise in the collisions. This scheme is characterized by a nonlinear 
numerical dissipation term with an automatic feedback mechanism which adjusts the amount of dissipation 
in any cell of the mesh according to the shape of the actual solution. The spatial discretization adopted in 
our code is based on upwind-differencing and the solution is advanced in time using a Godunov-type method 
(solving a set of Riemann problems at any cell interface forward in time). To obtain a 2" d -order spatial 



10 



accuracy we have used the TVD-MUSCL (Total Variation Diminishing-Monotone Upstream Scheme for 
Conservation Laws) reconstruction technique. We achieve 2 nd -order accuracy in time by using a two-step 
Runge-Kutta explicit scheme. 

The simulation starts when the two clouds are in contact. The cloud and intercloud density and 
temperature are computed from the corresponding phase diagram, as described in Sec. 2. The external 
medium is initially put into motion, with the same velocity as the clouds, to prevent the formation of a 
vacuum zone behind the cloud boundary due to the snowplow effect of the cloud motion. 

In Figure || we show, as an example, the evolution of the collision of two identical clouds for the 
representative case v r ~ 19 km s , L c = 0.1 pc, Z = D/G = 1. The three different phases discussed above 
are clearly seen in the various panels. At t — 0.5t co u two shocks are already well formed and propagate 
through the cloud compressing the gas. The pressure jump across the shock is oc v^ and the density jump 
is ~ 4 because the shock is still adiabatic. Then, the temperature of the shocked gas starts to decrease due 
to radiative losses occurring in the transition region (see Sec. 3). At t ~ t cM an expansion phase begins 
driven by the hot gas near the cloud boundary whereas the gas in the cloud interior is almost at rest. The 
shock propagates through the intercloud medium increasing both its pressure and temperature. The cloud 
expansion is finally halted by the enhanced external pressure. 

The behavior of the collision between two different clouds is initially analogous to that of identical 
clouds (except that the CD is not at the rest in the LSR). As the shock exits the smaller cloud, thus 
triggering its expansion phase, the ram pressure of the larger cloud accelerates the CD and the velocity of 
the shock in its interior decreases consequently 

Finally, we have performed some simulations of collisions between identical magnetized clouds with the 
field lines parallel to the shock front. As expected, the main effect is a lower compression of the postshock 
gas, which enhances the value of N coo i. This has important consequences on the elasticity of the collision as 
we will see in the next Section. 



5. Elasticity of the collisions 

In this Section we present (Figs. 3-7) the results concerning the elasticity of the collisions derived by 
analyzing the numerical simulations. The numerical results are compared with the appropriate analytical 
expressions obtained in Sec. 3. In order to achieve the best agreement between the two sets of results we 
have fixed the value of the free parameter a, introduced in eq. (H), to 1/3. With this choice we find a very 
good agreement also varying the cloud size and gas metallicity. The error on the numerical calculation of 
E r is less than 1% for all the runs, and it is due to the fact that we integrate the radiative losses during the 
simulation up to the time at which the fractional variation of the integral is < 1%. 

In the following we will discuss the detailed dependence of E r /E k on the various collision parameters, 
namely v r , L c , Z, k, /?, B. 

•E r /E k (v r ) (Fig. 3) If the relative velocity of the clouds, is v r -C v m or v r 3> v m , with v m ~ 30 km 
s , the collision is approximately inelastic {i.e. all the initial kinetic energy of the clouds is dissipated 
radiatively) . In a collision between identical clouds, this occurrence does not depend on any other parameter 
of the collision. For v r in the interval (v m — Sv,v m + Sv) where Sv ~ 10 km s , the energy dissipation 
in the collision decreases and for v ~ v m it reaches a minimum (i.e. the elasticity is maximum). The 
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value of v m is slightly dependent on the metallicity Z of the gas. This behavior of the elasticity as a 
function of the relative velocity can be understood as follows. We recall that N coo i = nov s t coo i, with 
v s ex v r ; since t coo i ex Vg/[noA(v s , Z)], then N coo i is the product of two terms: the first, uf, increasing, the 
second, A _1 (w s , Z), decreasing with v s . Thus N coo i has a maximum and the same is true for e ex N coo i. 
Differentiating epsilon with respect to v r we find that e is maximum when: 

9AK)_ dMT 1 _3_ 

dv r ~ 6 ' ° r 8T ~2' l bJ 

•E r /E l k (L c ) The elasticity of the collision increases when the cloud size is decreased because 7/ ex L^ 1 (Fig. 
3). 

•E r /E l k (Z) The elasticity of the collision increases when the gas metallicity is decreased because rj ex A -1 
and A is directly proportional to Z, as shown by Fig. 4. Interesting enough, collisions between clouds 
for which the relationship Zh\ rs const, holds, are equally elastic, for a given value of v r - The physical 
explanation can be found from an inspection of eq. ( |l0| ) and recalling that, to a first approximation, A ex Z 
and the equilibrium density no of the cloud is ex Z^ 1 ' 2 . This relationship is illustrated by Fig. 5. 
•E r /El(n, (3) When the cloud size ratio k = L\jL<2, and the cloud velocity ratio f3 = Vi/v% are not equal to 
one (i.e. different clouds), the maximum kinetic energy that can be dissipated in the collision is equal to 
the kinetic energy relative to the center of mass reference frame. Hence, the collision between two different 
clouds will be, in general, less dissipative with respect to the identical cloud case. Moreover, while the 
shock velocity in the smaller cloud is always ex v r , the shock will be slowed down in the larger cloud for the 
reasons discussed in Sec. 4. This will produce an additional minimum of E r /E % k at v > v m as shown by 
FigJq. The analytical approximation in that region tends to overestimate the elasticity by about 10 
•E r / E l k (B) The effect of an uniform magnetic field Bq in the cloud, is to decrease the energy dissipation 
(i.e. increase the elasticity) of the collision. Since the magnetic pressure limits the compression of the 
postshock gas; as a consequence, iV coo ; increases along with the elasticity of the collision (see Fig.0). In 
Figure we show the behavior of (1 — e) vs. B = Pb(0)/P(0) (the ratio of magnetic and thermal pressure 
at t = 0). 



6. Summary and discussion 



We have presented a first step towards the study of the collisions between diffuse interstellar clouds in 
a multiphase medium, with particular focus to the dissipation of the kinetic energy of the clouds (described 
through the elasticity parameter, i.e. the ratio of the final to the initial kinetic energy of the clouds) as a 
function of the parameters of the problem such as the cloud relative velocity, gas metallicity, cloud mass 
and velocity ratio, and magnetic field strength. 

The problem has been studied both analytically, obtaining approximate, albeit handy, relations valid 
for a wide range of parameter variation, and numerically, by means of a ID shock capturing TVD numerical 
code. The comparison between the two approaches has been explored in detail (see Figs. 3-7); we conclude 
that the agreement is very good. 

The following points summarize our main results: 
• The kinetic energy dissipation in cloud collisions is minimum (i.e. the collision elasticity is maximum) for 
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a cloud relative velocity v r ~ 30 km s -1 . 

• The above minimum value is proportional ZL\ (where Z is the gas metallicity and L c is the cloud size): 
the larger is Zl? c the more dissipative (i.e. more inelastic) the collision will be. 

• We find that the energy dissipation decreases (i.e. elasticity increases) when the magnetic field strength, 
and mass ratio of the clouds are increased and the metallicity is decreased, respectively. 

The previous results have been obtained assuming: (i) centered collisions: we have investigated a simple 
ID model of the collision only; the extension to the analogous 2D problem will be presented in a future 
paper. However, we have already explored the effects of the lateral outflow analytically, (ii) ionization 
equilibrium: this holds only approximately if the temperature of the post-shock gas if T < 3 x 10 4 K and 
Z > 0.01 (iii) the ionization of the pre-shock gas by the radiative precursor (Shull & McKee 1979) which 
could modify the shock structure is neglected, (iv) clouds are not self-gravitating. 

Our results, mainly aimed at deriving the dependence of the elasticity on a wide range of parameters 
characterizing the initial conditions of the collision in a multiphase medium, do not allow us to draw any 
firm conclusion on the final fate of the remnant of the collision (shattering or coalescence). In order to make 
progresses a 2D study is required. Nevertheless, our perception from the ID study presented here is that 
the presence of a magnetic field must be invoked in order to prevent the shattering of the clouds found by 
Klein et al. (1995). A related point is the possible phase transition from the CNM to the WNM associated 
with the collisions, which could result in a net mass and energy exchange from the cold to the warm phase. 
This could be relevant as far as the overall thermal equilibrium of the ISM is concerned. 

Cloud collisions might be responsible for the buildup of the observed mass spectrum (N(m) ex to~ 2,14 , 
Dickey & Garwood 1989) of diffuse clouds and for the formation of molecular complexes. Das & Jog (1996) 
study the evolution of molecular clouds under the effect of collisions and local gravitational interactions. 
Their recipe for the collision is a simple extrapolation of the results by Hausman (1981), i.e. coalescence 
occurs in a subsonic collision, fragmentation is the product of a supersonic one. However, most of the 
standard buildup models of diffuse clouds - inspired by the pioneering work of Field & Saslaw (1965) - 
assume that collisions are inelastic, a hypothesis that is too simplistic by far, as already made clear by our 
ID results. Also, as Jungwiert & Palous (1996) have pointed out, the elasticity of cloud collisions is a key 
parameter for the process of formation of multiple rings in disk galaxies. 

In a cosmological context, collisions between primeval clouds could be important for two reasons. On 
the one hand, the nonlinear evolution of primordial fluctuations is thought to generate a very clumpy state 
of the intergalactic medium. One can expect these protogalactic seeds to move through the background 
gas, collide and eventually coalesce; during this process part of their orbital momentum can be transformed 
into spin momentum of the merger (Chernin 1993). On the other hand, energy dissipation by means of 
cloud collisions in a forming galaxy is found to produce flattened systems (Gott & Thuan 1976, Larson 
1976, Abadi et al. 1990) and the disks of spiral galaxies. It is thus important to fully understand how the 
energy dissipation is affected both by the gas metallicity, which is supposed to be rather low, and by the 
phase structure of the ISM of these primeval objects. Our calculations provide a necessary ingredient for 
this type of studies. 



We are grateful to S. Balbus, J. Dickey, G. Field, R. Klein, C. McKee, and M. Shull for stimulating 
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Dominant cooling processes 


Notes 


Rcf. 


Cooling by fine-structure lines: 






CII(158^m) 


Impacts with H° and e~ 


1 


OI(63;44^m) 


Impacts with H°,e~ and H + 


2 


FeII(26;15/im),SiII(34.8Mm) 


Impacts with H° and e~ 


3 


Cooling by metastable lines: 


Impacts with e~ 


3 


CII(2326l) 






OI(6300;6363A) ,OII(3229;3726A) 


01 most important coolant 




NI(1.04 J um),NII(6548;6583A) 






FeII(5.34;4.12;1.26/«n) 






SII(6731;6717l) 






SiII(2240l) 






Cooling by highly ionized elements: 


Impacts with e~ 


4,5 


CIII,CIV,OIII,OIV,OV 


Important at T > 10 5 K 





Table 1: Dominant cooling processes by metal lines. Reference: (f)Wolfire et al., 1995;(2)Pcquignot, 
1990;(3)Hollenbach & McKee, 1989;(4)Gaetz & Salpeter, I983;(5)Shull & van Steenberg, 1982. 
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Fig. 1. — Thermal pressure P/k vs. hydrogen density n; the curves refer to different values of the dust- 
to-gas ratio, D/G, and metallicity Z, with D/G = Z. The gas is thermally stable for (d\ogP/d\ogn) > 
(i.e. positive slope of the curves). Unless Z = 0, a stable two-phase medium is supported. 



Fig. 2. — Evolution of a collision between identical clouds with relative velocity v r — 19 km s -1 , size L c = 0.1 
pc, Z = D/G = 1. The panels show the density, velocity, pressure and temperature at the evolutionary 
times t ~ 0.5,0.98,1.6,4 t colli with t co ll = L c /v c ~ 10 4 years. 



Fig. 3. — Behavior of the quantity (1 — e) = E r /E k , where e is the elasticity of the collision (see text) for 
a collision between two identical clouds of size L c shown by the labels, as a function of the cloud relative 
velocity v r . The gas metallicity is Z = 1. The solid curves represent the analytical expression for 1 — e 
(where a = 1/3; see text); the points show the analogous results from the numerical simulations. Vertical 
bars denote the estimated error in the simulation (see text). 



Fig. 4. — Same as Fig. 4, but for different values of D/G, and gas metallicity Z. We assume D/G = Z, 
L c = 1 pc and Tq = 48 K; the cloud density no is derived from the corresponding phase diagram. Vertical 
bars denote the estimated error in the simulation (see text). 



Fig. 5. — Contour levels of the function 1 — e = E r /E l k in the plane L c — Z for a relative velocity v r = 30 



Fig. 6. — Same as Fig. 4, but for a collision between different clouds as a function of the cloud relative 
velocity v r ; the curves refer to different values of the cloud size K — L c \/ L C 2 and velocity (3 — | Vi /^2 1 ratios; 
we have assumed (3 2 — 1/k so that the clouds have the same initial kinetic energy. The points show the 
analogous results from the numerical simulations. Vertical bars denote the estimated error in the simulation 
(see text). 



Fig. 7. — 1 — e = E r /E k curves plotted as a function of the initial (t — 0) ratio between the magnetic 
and thermal pressure Pb(0)/P(0) in the cloud. The points shown in the picture are computed for: 
Pb(0)/P(0) = 0.1; 0.5; 1; 5; 10 (i.e. B ~ 1; 2; 3; 6; 9^G). The collision parameters are the same as in Fig. g. 
Vertical bars denote the estimated error in the simulation (see text). 
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